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introduction 


A.  Integration  of  Single  Rays 

We  are  dealing  with  the  ray  approach  to  mapping  sound  propagation 
in  the  sea . 

We  are  considering  the  problem  in  which  the  speed  of  sound,  c , is 
variable  over  two  dimensions,  the  depth,  z , and  a horizontal  coordinate,  x . 
With  no  refraction  normal  to  the  x,z  plane,  a ray  directed  in  the  plane 
remains  in  the  plane. 


0 x 


fig.  1 A ray  is  the  path  of  an  element  of  wavefront 


Each  element  of  wavefront  is  instantaneously  propagating  along  its 
normal,  in  the  positive  direction  denoted  by  unit  vector  It  . A ray  may  be 
defined  as  the  path,  or  locus,  of  an  element  of  wavefront.  The  path  is 
bent  by  refraction.  The  governing  equation  may  be  expressed  by 


& 

ds 


COS  > 

R - z 


i *£ 

c ?»n 


(1) 


The  notation  is  identified  with  reference  to  Fig.  1:  The  direction  It  is 
represented  by  the  angle  y measured  clockwise  positive  from  the  horizontal. 
Linear  distance  measure  along  the  ray  is  denoted  by  s . Linear  measure 
normal  to  the  ray,  and  hence  tangent  to  the  wavefront,  along  unit  vector  n 
is  denoted  by  n . The  unit  vector  n is  clockwise  normal  to  It  by  sign 
convention. 

The  first  term  on  the  right  hand  side  of  Eq.  (1)  is  due  to  coordinate 
curvature  of  the  reference  horizontal;  R is  the  radius  of  curvature  of  local 
sea  level.  The  second  term  expresses  the  refraction,  caused  by  speed 
gradient  along  the  wavefront. 

R . Elimination  of  the  Curvature  Term 


Equation  (1)  may  be  transformed  into 


d y 

ds 


Sn 


U) 


where 


c 


* 


(3) 


and  x is  an  arbitrary  constant.  We  choose  to  make  x = R. 

The  factor  R/(R-z)  ranges  from  1 at  the  surface  to  1.001  at  a depth 
of  about  6 km.  Over  this  range  of  depth  the  factor  increases  the  sound 
speed  almost  linearly  with  depth  by  zero  at  the  surface  to  about  1.5  meters 
per  second  at  6 km.  The  significance  for  including  the  curvature  effect 
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may  bo  judged  In  the  context  of  the  variability  of  the  sound-speed  structure 
and  the  accuracy  with  which  it  may  be  specified  from  available  information. 
In  any  case  the  curvature  is  a bias  which  may  readily  be  Included  by 
modifying  the  speed  according  to  L’q.  (3). 

The  linear  argument,  s , along  a ray  may  be  replaced  by  travel 
time,  t , according  to 


ds  = c dt 


The  use  of  c*  in  place  of  c results  in  a pseudo  travel  time 


= c*  dt* 


dt 


* 


dt 


(l) 

(5) 

(6) 


The  modification  is  slight. 

In  all  that  follows  we  shall  consistently  use  the  modified  speed, 
and,  in  the  few  places  it  appears,  the  modified  travel  time.  Hence  we 
can  omit  the  subscript  asterisk  without  confusion.  The  ray  tracing  is 
governed  by 


±2  x _ i (7) 

ds  c 3n 

C.  The  Specific  Wavefront  Length 

An  individual  ray  may  be  identified  at  a reference  point  denoted  by 
s in  Tnj.  I,  where  the  orientation  of  the  particular  wavefront  element  is 
defined  by  V , . The  wavefront  element  also  has  two  dimensions  associated 

s 


- « 


l' 


-3- 


with  its  stretching  (wavefront  divergence)  or  shrinking  (wavefront  conver- 
gence) and  reversal  (wavefront  folding) . We  denote  the  specific  wavefront 
length  in  the  x,z  plane  by  L,  and  shall  not  concern  ourselves,  for  present 
purposes,  with  the  lateral  dimension  of  the  wavefront. 

The  derivation  of  the  governing  equation  for  L along  the  ray  is  given 
in  Reference  [l]  . This  governing  equation  may  be  expressed  by 

d^L  1 dc  dL  L d2c 

2 c Ss  ds  c 2 

ds  an 

The  specific  wavefront  length  may  be  initialized  (normalized)  at  the  point  S 
by  specification  of 


Initialization  also  requires  specification  of  the  curvature  of  the  wavefront 
at  S , by  specification  of 


(dL/ds)g 

Equation  (8)  is  for  single-ray  integration  of  wavefront  spreading. 

D.  Depiction  of  Wavefront  Propagation 

The  ray  approach  may  be  used  to  depict  the  propagation  of  wave- 
fronts  by  integration  of  the  traces  of  sufficient  rays  to  adequately  depict 
the  wavefront  continuum  at  all  positions  in  its  propagation.  Such  complete 
mapping  of  the  wavefront  propagation  encompasses  the  wavefront  spreading, 
shrinking  and  foldings,  and  implicitly  includes  full  determination  of  L. 


Independent  integration,  along  some  ray,  of  L by  a numerical 
analogue  of  Eq.  (8)  thus  yields  another  determination  of  L.  Because  of 
differences  in  truncation  errors  between  the  two  independent  numerical 
integrations--rays  by  numerical  analogue  of  Eq.  (7)  and  L by  numerical 
analogue  of  Eq.  (8) — the  results  may  not  be  consistent.  For  complex 
sound  distributions  or  badly  designed  numerical  analogues  the  results  may 
not  even  match  in  areas  of  interest.  In  any  event,  because  truncation 
errors  compound  with  length  of  the  independent  integrations,  the  match 
deteriorates  with  distance  covered  by  the  rays. 

Of  pragmatic  interest  is  the  propagation  from  a point  source.  Let  S 
of  Fig.  1 be  the  location  of  the  point  source.  Then  each  ray  is  defined  by 
its  emission  angle,  yg  , and,  for  all  rays, 

Ls  , 0 01 

(dL/ds)s  = 1 . (10) 

Any  point  in  the  medium  may  be  traversed  by  one  or  more  rays, 

giving,  in  association,  to  that  point  one  or  more  values  of  yg  . Such 

multivaluedness  is  produced  by  foldings  of  the  wavefront.  The  wavefront 

remains  continuous  if  we  include  the  foldings  in  our  concept  of  continuity, 

and,  correspondingly,  the  field  of  >g  , while  multivalued,  is  also 

continuous.  Thus  multivalued  areas  are  separable  into  superimposed 

families  of  continuous  y distributions.  For  each  family  the  ascendent 

o 

of  y is  tangent  to  the  wavefront  in  the  direction  of  increasing  source 
S 

emission  angle.  This  direction  is  + n or  - in  depending  on  the  even- 
or-oddness  of  the  number  of  foldings  the  wavefront  element  has  undergone 
in  its  propagation  from  the  source. 
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For  a point  source  the  specific  wavefront  length,  L , may  be 
diagnosed  from  full  ray  depiction,  throughout  the  insonified  medium,  by 

f = V (11) 

for  each  family.  Accordingly,  and  alternately  but  consistent  with  our 
earlier  definition,  L may  be  defined  as  the  specific  wavefront  length  per 
unit  radian  of  emission  at  tire  point  source.  This  definition  is  restricted 
to  insonificatlon  by  a point  source,  whereas  the  earlier  definition  is 
general . 

F.  Preference  for  Consistency 

In  mapping  the  sound  propagated  from  a point  source  we  have 
available  two  methods  for  determining  the  multi-valued  field  of  L.  The 
first  method  involves  tire  design  and  application  of  a numerical  analogue 
for  the  ray  tracing  equation,  Eq.  (7),  and  for  diagnosis  of  the  ray  spacing 
by  Eq.  (11).  The  second  method  involves  the  design  and  application  of  a 
numerical  analogue  for  Eq.  (it)  for  integration  along  each  ray  obtained  by 
analogue  of  Eq.  (7).  Unless  the  numerical  analogue  for  Eq.  (8)  can  be, 
and  is,  designed  in  a form  consistent  with  the  ray  tracing,  truncation 
errors  will  manifest  quite  differently.  Ihe  disparity  will  generally  grow 
with  range  and  with  the  complexity  of  the  sound-speed  distribution  and 
the  ocean-bottom  topography. 

Such  independence  in  truncation  error  may  be  deemed  useful  as 
indicator  of  the  range  limits  imposed  by  truncation  errors  overwhelming 
resolution.  However  we  dismiss  this  use  because  the  likelihood  is  that 
the  ray  spacing  would  be  more  accurate,  in  general. 


We  deem  it  dc  sirable  to  compute  specific  wavefront  length.  L , 
rather  than  to  rely  on  spacing  between  rays  of  finite  separation  and  diverse 
paths.  And  we  deem  it  desirable  that  L be  consistent  with  the  distribution 
of  rays . 

Wc  consider  the  information  to  be  complementary.  Continuity  in 
L from  ray  to  ray  is  an  indicator  of  the  adequacy  of  the  resolution. 

The  ray  to  ray  locus  of  L - 0 locates  caustics.  Caustics  generally 
do  not  lie  at  the  intersection  point  of  two  rays  of  finite  separation  in 
but  occur  at  the  differential  crossing  or  folding.  Propagation  loss  by 
diffraction  depends  on  the  geometry  of  the  propagation  with  strong 
dependence  on  the  occurance  and  orientation  of  caustics. 

F.  Objective 

Our  objective  is  to  develop  a numerical  scheme  for  ray-consistent 
determination  of  the  specific  wavefront  length  distribution.  This  distribution 
generally  becomes  more  multivalued  with  range  from  the  source,  not  only 
by  internal  foldings  in  sound  channels  but  also  by  reflections  at  sea  surface 
and  bottom. 

This  endeavour  involves  considerations  which  include 

— maximizing  accuracy  in  ray  tracing  and  ray  spacing  for  the  prescribed 
sound  speed  distribution, 

— simplifying  the  integrations  for  the  purpose  of  reducing  calculations 
for  the  same  resolution  yield, 

— exploiting  simplifications  which  occur  when  sound  speed,  c , is  a 
function  only  of  depth,  z , and 

— realizing  consistency  as  horizontal  variability  in  c vanishes. 

We  proceed  by  expanding  on  aspects  of  single-ray  integrations  and 
simplifying  the  equations  for  when  speed  is  a function  of  depth  only.  We 
then  return  to  designing  the  object  scheme. 
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II 


Single  Ray  Spreading  Integration 


A.  Transformations  of  the  Governing  Equation 

Equation  (8)  is  a form  of  the  differential  relationship  which  governs 
the  specific  wavefront  length,  L , along  a ray.  The  motivation  for  trans- 
forming this  equation  into  other  forms  includes  revealing  the  character  of 
the  equation,  the  significances  in  speed-of-sound  distributions,  and 
insights  for  designing  improved  numerical  analogues. 

Special  treatment  of  the  equation  is  required  for  use  in  integrating 
through  discontinuities  in  speed  distribution  at  interfaces  and  at  surface 
and  bottom  reflections.  This  development  follows  in  Subsection  B. 

We  begin  the  transformations  by  introducing  K 5 C*:c  . 


d2L  1 dL  ( S2K 

. 2 ~ 3s  ds  j ^ 2 

ds  l 3n 

Keep  in  mind  that  s is  along  the  curved  path  of  the  ray  and  n is  along  a 
straight  line  normal  to  the  ray.  At  a point  along  the  ray  we  can  write 

V2K  = 32K/3n2  + 32K/3s*2  (13) 


where  s^  is  along  a straight  line  tangent  to  the  ray.  We  also  require  the 
operator  transformation 


a 2 

■—  s it  • v (it  • V)  = (tt  • v)  + tt  • v tt  - v 

3s2 

32  dy  32  _ 3K  _3_ 

^2  111  ds  ,2  3n  3n 

as*  Ss* 


(14) 
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In  the  last  step  we  have  used  Kq.  (7).  Accordingly  wo  have  the  trans- 
foi  mation 


► 


A*\k 


(15) 


By  using  Lqs.  (H)  and  (15)  we  tiansform  Kq.  (12)  into 


d l dK  dL 


l’his  tor m is  significant  because  the  directional  operator,  ( *n  • ^ ) , ol 

> 

Eq.  (B)  has  been  leplaced  in  favor  ot  an  isotropic  operator,  v“,  and 
differentials  along  the  ray. 

In  Reference  [2]  we  have  transformed  the  governing  equation  into 
the  toi  in 


I v2 

\ v K 


A { 


(16) 


) > 

- c“  V“K  Q 


(17) 


where  0 1.  c and  travel  time,  t , replaces  s as  measure  along  the  ray 

according  to  ds  cdt.  This  is  the  simplest  fotin  we  have  been  able  to 
deduce.  It  enabled  us  in  Reference  [2]  to  determine  that  0 is  reciprocal 
between  a point  source  and  a receiver,  for  any  ray  path. 

Initializations , at  a point  source,  are  given  by  Eqs.  (9)  and  (10) 
and  theii  transforms: 

Qs  0 ; (d Q dsK  1/c  ; (dQ  dt)^  1 (1^) 
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B.  Special  Treatment  at  Discontinuities  and  Reflections 


In  describing  the  speed  distribution  in  the  sea  it  nuiy  be  convenient 
or  suitable  tv'  include  interfaces  where  the  speed  gradient  abruptly  changes 
rot  present  purposes  we  shall  treat  such  interlaces  as  horizontal.  Ray 
tracing  by  Eq.  (7)  requires  only  that  the  speed,  c , be  single-valued  and 
continuous  (zero-order  continuity).  It  can  be  integrated  through  interfaces, 
of  abrupt  change  in  speed  gradient  (first-ordei  discontinuity) . llowevei 
Eq.  (B)  cannot  be  integrated  through  such  interfaces;  the  spreading  rati', 
d I ds,  abruptly  changes . 


Reflections  at  the  sea  surface  or  ocean  bottom,  where,  m general, 
we  wish  to  allow  3c  3.-  * 0 , also  involve  effective  passage  through  kink: 
in  the  speed  profile.  At  reflections  we  have  the  additional  consideration 
that  1 , in  accordance  with  Eq.  (11),  tlips  sign. 

The  expressions  required  for  carrying  the  integration  of  Eq.  (8) 
through  interfaces  and  reflections  have  been  derived  m Reference  [ 3 . 

This  was  done  by  applying  Eq.  (8)  to  a transition  layer  between  levels  ot 
differing  specified  gradient;  letting  the  thickness  ot  the  layei  shrink  to 
zero  in  the  analysis  gave  the  desired  expression. 


At  reflections  the  sign  of  1.  is  flipped. 


(It! 


1. 


R 


cc 


The  subscript  1 is  used  tv'  denote  incident  value  and  the  subscript  R a 
reflected  value.  The  spreading  rate  changes  according  to 


Sc 

Sz 


(21) 


2 

2Lr  cos  Vj 
c sin 


Transformation  of  these  equations  in  terms  of  Q is  readily  effected 
by  substitution  of 


L = Qc  ; dL/ds  = c dQ/ds  + Q Sc/Ss 


(22) 


Sc/Ss  = cos  y dc/Sx  + sin  y Sc/Sz 


(23) 


The  transformation  yields 


C.  Numerical  Analogue  Design  and  Accuracy  Demands 

Consider  the  design  of  a numerical  analogue  for  the  spreading 
integration  along  a single  ray,  without  any  regard  for  consistency  with 
ray  spacing.  The  simplest  design  is  suggested  by  the  Eq.  (17)  form  of 
the  governing  differential  equation.  This  design  is 
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Section  II  B that  sharp  changes  in  profile  gradient  cause  pronounced  changes 
in  the  spreading  rate.  Unless  such  features  are  treated  as  kinks,  Eq.  (~B) 
must  be  Integrated  in  very  small  increments,  through  these  features,  In 
order  to  contain  truncation  errors. 

This  sensitivity  would  be  improved  if  the  numerical  analogue  could 
be  expressed  in  terms  of  perfect  differentials,  at  least  in  the  depth 
component  ot  the  argument.  Reduction  ol  the  differential  order  oi  the  speed 
dependence  acts  m the  same  direction  of  reducing  sensitivity. 

Either  success  would  also  overcome  the  need  for  special  treatment 
at  internal  profile  kinks  and  at  reflections  (provided  we  flip  the  sign  ot  1.). 
The  consistent  scheme  that  we  outline  in  Section  V has  such  advantages. 
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1 1 1 Snujlt'-rioHU*  Stmpli  beat  tons 

A.  Ray  Tracing 

lu  the  present  section  we  treat  the  special  case  of  negligible 
horizontal  variation  in  sound  speed.  The  speed,  c , is  prescribed  as  a 
function  only  ot  depth,  z . 

liquation  (7),  for  a ray,  may  be  written 

dT  _ cos  "Y  dc  (11) 

ds  = " c dz 

Snell's  law  may  be  obtained  by  considering  differentials  for  an  increment  in 
depth,  flz  . 


ft) 


- cos  y fti  ■ c 


6s 

6z 


By  substituting  sin  ) ft 6s  we  obtain 

6 . ■ cos  y ft  «.’■  c 


which  means  that 


Is  constant  along  a ray. 


(d  J) 


Retelling  its  value  to  the  leference  point  S we  may  write, 


This  is  Snell's  law,  and  we  have  shown  that  it  can  include  the  horizontal 
coordinate  curvature  factor  by  modification  of  the  speed  profile. 


Snell's  law  may  be  useful  in  ray  tracing.  It  fails  to  span  extremal 
depths  (i.e.  over  > = 0).  This  limitation  arises  because  multiplication  by 
sin  y is  implicit  in  its  derivation  from  Eq.  (31) . 

Let  us  denote  successive  values  along  a ray  by  subscripts  .... 
m-1,  m,  m+1  . . . Equation  (32)  yields 


cos 


m+1 


Cm+ 1 
c 

m 


cos  y 

m 


(34) 


A form  of  centered  scheme  is  expressed  by 


6x  = 4 (cos  y + cos  y , ) 6s 

Z m m+ 1 

6z  = 4 (sin  y + sin  V , ) 6s 

2 m m+ 1 


(35) 


(36) 


If  we  choose  6s  as  increment  then  Eqs.  (34),  (35)  and  (36)  are  an  implicit 
set  requiring  iteration.  If  we  choose  6z  as  increment  their  use  is 
explicit. 

With  6z  as  increment  the  integration  proceeds  nicely  until 


cos  y , > 1 

m+1 

For  such  increments,  m to  m+1,  we  must  allow  6z  to  be  determined  by 

cos  y = 1,  and  the  progression  in  z is  reversed  for  the  subsequent 
m+ 1 

step.  At  reflections  we  must  also  allow  6z  to  be  suitably  specified. 


B.  The  Spreading  Equation 


We  specify  a point  source.  Snell's  law,  expressed  by  Eq.  (33), 
holds  everywhere.  The  specific  wavefront  length  and  the  spreading  rate 
at  a point  on  any  ray  are  given  by 


i 

L 


(37) 


i dL  = ay 

L ds  9n 


(38) 


We  perform  the  operation  S/Sn  on  Snell's  law  at  the  arbitrary  point,  and 
obtain 


sin  y by 
c bn 


cos  y Sc 
2 Sn 


sin  ST, 


Sn 


(39) 


Substitution  according  to  Eqs.  (37)  and  (38)  leads  to 


dL 

ds 


sin  ys 

sin  y 


c 

Co 


L 

c 


2 V 

cos  y 
sin  y 


Sc 

Sz 


(40) 


This  is  an  order  lower  than  the  general  form,  Eq.  (8).  Equation  (40)  fails 
to  integrate  through  points  at  which  y = 0 , but  gives  the  result  at  these 
points  that 


2 sin  Tc 
T = c_  S 

Cg  Sc/Sz 

This  suggests  that  Eq.  (40)  may  have  more  to  offer. 


(41) 
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Equation  (40)  may  be  cast  Into  a simpler  form  by  substitution  of 


L = X sin  y (42) 

The  new  principal  dependent  variable,  X , may  be  interpreted  as  the 
horizontal  component  of  the  ray  spacing.  Substitution  leads  to 


dX 

ds 


tan  >s 


cos  y 
2 

sin  y 


which  may  also  be  written 


dX 

ds 


cssinys  c 

2 2 2 
cg  - c cos 


(4  3) 


(44) 


This  reveals  that  dX/ds,  along  a ray,  is  a function  of  c(z)  only.  We  shall 
use  this  fact  to  prove  recursion  formulae  along  rays.  We  shall  use  the 
schematic  form  of  Eq.  (44)  , 


^ = F(z)  (45) 

ds 


C.  Recursion  Formulae 

We  are  still  dealing  with  a single-profile  speed  distribution,  c(z). 
For  a ray,  defined  by  y g , Eq.  (45)  may  be  written 


dX  = __F 

dz  sin  y 


(46) 
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(4  8) 


2\ 


In  segment  3 we  have 

X, 


*b  + 


.1  sin 


y „ 


dz 


*b  + 


J sin 


dz 


In  segment  4 we  have 


X = \ + ^ 

4 c J 


sin  y 


dz 


4 


= \ + \ y 

c j Sin  y 


dz 


- S 


sin  y , 


dz 


= X + 
c 


sin  V, 


dz 


- s 


sin  y„ 


dz 


- 2X 


X . = 2X 

d c 


(49) 

(50) 


Subsequent  cycles  have  the  additive  values:  X^  , 2X^  , 3X^  , etc. 

These  formulae  were  noted  and  exploited  by  E.  Hesse  of  FNWC. 
The  present  verification  followed. 
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IV 


A Scheme  for  Ray  Tracing 


Tor  the  general  case,  in  which  c varies  over  the  x,z  plane,  we 
would  like  a scheme  that  affords  consistency  with  Snell's  law  as 
dc/dx  ■>  0. 


We  are  inclined  to  prefer  6z  as  the  independent  increment  for  the 

ray  integration,  except  for  turning  levels  (V  0)  where  6z  is  determined 

to  the  y - 0 level.  Also  the  increment  can  be  adjusted  to  accommodate 

the  z levels  of  interfaces  and  surface  and  bottom  reflections.  Another 

motivation  in  specifying  6z  is  that  this  increment  can  readily  be  geared 

to  the  variability  of  the  speed  structure.  For  example  consider  making  6z 

2 . 2 

inversely  proportional  to  d c/dz  . 


By  noting  that 


Ss 


cos  y 


3x 


+ 


sin  y 


dz 


(51) 


_a_ 

3n 


- sin  y T- 
?x 


cos  y 


?z 


(52) 


and  K c , we  transform  Eg.  (7)  as  follows: 


ill 

ds 


sin  y 


dK 

c>x 


cos  y 


?z 


sin>  dy  l-cos°y  dK  . v dK 

— = r — ~ + sin  7 — 

cosy  ds  cosy  ox  oz 


_d_ 

ds 


ir.  cos  y 


1 jVK 
cos  y dx 


d cos  y _ 1 

ds  k c cos  y dx 


(53) 
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- -1 


. 


1 


4 


We  may  develop  the  numerical  analogue  in  terms  of  successive 
points  along  the  ray,  . . . . m-1,  m,  m+1,  . . . Exact  application 
yields 


m+  1 

-5 


1 9K 

t — ds 


cos  * 


m+  1 


m+1  ..  m 

cos  y e 

c m 

m 


cos  y ds 


(54) 


which  we  approximate  by 


cos 


s y 


m+  1 


m+  1 


m 


cos 


T,  6s  3K 

y 1 - v T- 

m cos  y ds 


(55) 


It  remains  to  decide  how  to  evaluate  the  term  under  the  curly  line.  Since 
it  is  generally  acknowledged  that  the  term  is  very  much  smaller  than  1 
a forward  evaluation  should  be  adequate.  We  substitute 

6s  = 6z/sin  y 


and  obtain 


This  analogue  is,  however,  not  suitable  at  turning  levels  O'  =0).  We 
design  a special  operation  for  the  turning  increment. 


V 


■ 


1 


Fig.  3 Integration  through  V = 0 


For  increments  at  the  turning  level,  as  shown  in  Fig.  3,  we 
propose  a numerical  analogue  obtained  from 


d>  _ _ 3K  _ 1 

ds  dz  r 


(57) 


That  is,  we  approximate  the  trace  by  circle  arcs  with  radius  of  curvature  r, 
and  obtain 


♦ 

0 


(59) 


(60) 


ym+2 


(61) 
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4* 

T 


It  may  also  bo  strongly  argued  that  the  relative  contribution  of 
3K/3x  in  Eq.  (56)  is  negligible.  If  the  purist  does  not  agree  then  he  must 
admit  to  ranges  for  which  lateral  refraction  is  of  comparable  significance. 
At  large  ranges,  it  should  be  realized,  concepts  of  apertures  and  ducting 
replace  dependence  on  exact  ray  positions. 

In  any  case,  we  have  realized  consistency  with  Snell's  law  by  our 
numerical  analogue,  Eq.  (56).  The  ray  tracing  may  be  performed  as 
follows: 

(1)  Specify  fiz  for  the  increment  m to  m+1,  interpolate  c ^ and 
(3K/?x) 

m 

(2)  Compute  cos  > , from  Eq.  (56)  . 

m+1 

(3)  If  cos  > > 1 abandon  the  value  and  use  the  turning  operation 

m+1  — 

defined  by  Eqs.  (58)  through  (61)  to  perform  two  increments  to 
point  m+2  . 

(4)  If  cos  y < 1 , compute  sin  > . from 

m+ 1 m+ i 


1/2 


sin 


m+ 1 


= ( 1 - cos' 


W 


(62) 


with  positive  sign  for  downward  progress,  6z>  0 , and  negative 
for  upward  progress,  6z  < 0 . 

(5)  Compute 


fix 


6z 


cos  y + 
m 

sin  y + 
m 


cos  y 
sin  y 


m+  1 
m+ 1 


(63) 


The  angle,  y , need  not  be  evaluated. 

The  remaining  details  of  the  ray  tracing  integration  should  be 
straight  forward. 
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V 


A.  General  Concept 

The  governing  equation,  Eq.  (8),  for  the  wavefront  spreading  was 
derived,  in  Reference  [l]  , by  a combination  of  the  refraction  equation, 


d7  = _ 1 Sc 

ds  c 8 n 


and  the  spreading  rate, 


i dL  = 

L ds  9n 


(64) 


(65) 


The  propagation  of  a finite  segment  of  wavefront  can  be  depicted  by  an 
adequate  density  of  rays  each  integrated  by  Eq.  (64) . The  adjunct 
integration  of  the  spreading  must,  for  consistency,  be  performed  in  complete 
accord. 

We  propose  that  the  adjunct  spreading  integration  be  performed  in 
terms  of  ambient  segments  of  rays  for  each  increment  of  the  ray  tracing. 

The  scheme  is  depicted  in  Fig.  4. 


Fig.  4 Adjunct  Ray  Segments  Along  a Ray  Tracing. 
Normal  Segmentation 
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For  each  increment  of  the  ray  tracing  the  ambient  ray  segment  is 
initialized  at  a point  one  suitable  length  unit  removed  along  the  normal, 
on  the  appropriate  side.  The  side  depends  on  the  sign  of  L • 

The  Initialization  at  adjunct  point  m , where  Lm  and  (dl/ds)m 
are  known,  is 


y * = y + L 1 (dL/ds)  (66) 

m m m m 


The  ray  segment  is  calculated  by  the  identical  scheme  used  for  the  lay. 

Consider  extending  the  segment  to  the  next  ray  normal  at  point  m+1 ; it 

meets  this  normal  at  distance  6nm+1  from  the  point,  with  new  angle 

y*  . i'rom  these  values  compute 
m+  1 


L 


m+  1 


fin 


m+ 1 


L 

m 


(67) 


(db/ds) 


m+  1 


= L 


m+1 


m + 1 


rtvt- 1 


fin 


m+  1 


(68) 


Note  that  L may  change  sign;  the  subsequent  segment  then  switches 
sides . 

The  scheme  may  be  considered  a refinement  of  using  pairs  of  tays. 
It  is  an  improvement  on  pairs  because  pairs  may  diverge  considerably  and 
become  non-representative. 

This  scheme  can  be  integrated  through  kinks  in  the  speed  profile 
and  at  reflections , without  special  treatment  except  for  flipping  the  sign 
of  L at  reflections. 
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B.  Particular  Design 


Particular  design  of  the  adjunct  spreading  integration  must  be  tied 
in  with  the  scheme  adapted  for  the  ray  tracing.  We  will  now  go  along  with 
the  scheme  outlined  in  Section  IV. 

Depth  increment,  6z  , is  used  as  the  independent  increment  for  the 
ray  integration,  except  for  turning  levels  (V  = 0).  The  scheme  is  depicted  in 
Fig.  5. 


Fig.  5 Adjunct  Ray  Segments  Along  a Ray  Tracing. 

Horizontal  Segmentation 

As  principal  dependent  parameter  in  place  of  the  specific  wave- 
front  length  we  will  use 

X s L/sin  y (69) 

We  require  one  more  transformation.  Substitution  of  Eq.  (69)  into 
Eq.  (65)  yields 


i dA 

\ ds 


cos  v yy_ 
sin  y ds 


yy 

On 


;iu  > d\ 
\ cts 


"OS  ^ 


oy 

Os 


t sin  y 


a y 

dn 


Expansion  of  operators  by  Eqs.  (51)  and  (52)  leads  lo 


sin  y d\ 
\ ds 


Ox 


(70) 


A more  convenient  form  is 


0 cos  y 
ax 


sin  y d\ 
X ds 


(71) 


Each  adjunct  segment  is  performed  as  follows.  In  lie.  5,  at  point 
m-2  for  example  let  \ m_  , and  (d\/ds)m_2  already  have  been  computed. 
The  adjunct  segment  is  placed  one  suitable  length  unit  removed,  along 
minus  x toi  \ positive.  The  segment  is  directed  according  to 


'os  y 


cos  y 0 t 
m-2 


sin  > 

m-2 


dV 

ds 


(72) 


m 


_2 


The  segment  is  extended  over  Increment  6s*.  , the  same  increment  used 
foi  the  ray.  Equation  (5t>) , applied,  yields 


cos  y 


m- 1 


From  those, 


r 


I 

t 


\ i " A x , X •> 
ni-1  m- 1 m- 1 


(74) 


(dXAls)  - — j“‘C|-co,,..i|  (7i) 

sin  r . 
m-  1 


At  tiio  turning  levels  (>  0)  tin'  scheme  ditlers  in  that  tin'  adjunct  segment 

is  taken  ovei  two  lncreinents--l.e.  around  the  oxtremal--accordiny  to 
Fqs.  (S7)  through  (i>  1)  ns  used  tot  tlio  ray  itst'lt. 

Wt'  hnvt'  probably  overlooked  sonu'  details  .uni  complications  iiiny’ 
remain.  Howevci  we  tiust  that  they  can  be  ironed  out  in  programming  and 
testing . 
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